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Abstract 

By using geometric methods and superenergy tensors, we find new simple criteria 
for the causal propagation of physical fields in spacetimes of any dimension. The 
method can be applied easily to many different theories and to arbitrary fields (such 
as scalar or electromagnetic ones). In particular, it provides a conservation theorem 
of the free gravitational field in all ^-dimensional spacetimes conformally related to 
Einstein spaces (including the Ricci-flat case) . In general relativity, our criteria provide 
simple proofs and a unified treatment of conservation theorems for neutrinos, photons, 
electrons and all other massless and massive free spin n/2 fields. The uniqueness 
of the solution to the field equations also follows from our treatment under certain 
circumstances. 
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The traditional way of studying general properties of field equations for physi- 
cal fields on spacetimes is to introduce some coordinates and translate the geometric 
structure into the language of standard theory of partial differential equations to apply 
analytic methods or to consider the existence of fundamental solutions (propagators), 
sec [||. ||, H, [|, ||, ||, ^, P and references therein. In this letter we put forward a way to 
use the geometric information of tensorial field equations in a more direct manner. We 
show how a conservation theorem for fields can be proven, a theorem which provides 
information about the causal propagation of the fields. Our methods work for space- 
times of any dimension, so that it can be applied to many different theories including 
string theory, M-theory (see e.g. Q), supergravity, gauge theories, Kaluza-Klein-like 
theories, Brans-Dicke and scalar-tensor theories, dilaton gravity and many others. The 
idea is to generalize the conservation theorem proved in [10] for the energy- momentum 



tensors that satisfy the dominant energy condition to the case of arbitrary physical 
fields. The key points in our reasoning are: first, "superenergy" tensors with the cor- 
responding dominant property can always be constructed and, second, these tensors 
do not have to be related to any physically meaningful quantity, nor to the energy- 
momentum tensor of the field (which, in fact, cannot be used in many cases, see below, 
including the important case of the gravitational field where it cannot be defined) . The 
only condition in our conservation theorem is an assumption on the divergence of the 
superenergy tensor which is easy to check explicitly. 

Our method is clearly related to the techniques of symmetric hyperbolic systems. 



Friedrich j|, 11 1 has used these techniques to study many interesting equations including 
several of the examples we will discuss here. The great advantage of our method is 
that we do not need to put the equations on a form which fits the theory of symmetric 
hyperbolic systems but we use the superenergy tensor directly in a typical integration 
procedure and to obtain estimates. 

Let spacetime Vjv be an TV-dimensional Lorentzian manifold with a metric g of 
signature (+, -, -). A rank-m tensor T is said to satisfy the dominant property if 



El 0,0 



1 a 1 ...a m U\ ■ ■ ■ U m ^ u 



for any set {u\, . . . ,%} of causal future-pointing vectors. This is equivalent [14, 12] 
to that with respect to an orthonormal basis {e*o, . . . , e*7v-i} the pure time component 
dominates any other component, that is 



! o...o 



> \T a ...p\ (1) 



for any < a, . . . , (3 < N — 1. Note that this also implies that T0...0 = => T = 
0. (Observe that the dominant property is a generalization of the dominant energy 
condition for rank-2 energy-momentum tensors pi.) 
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In general, a physical field A will not have the dominant property but, as has been 
proven in complete generality fl^, [DJ and in arbitrary dimensions [|l4|] , one can always 
construct superenergy tensors T(A) of A which have the dominant property and such 
that, in any open neighbourhood, 

A = <=^T(A) = 0. (2) 

This property holds in open neighbourhoods of the manifold, but it may need a refine- 
ment at isolated points or lower-dimensional subsets. For instance, sometimes T(A) 
vanishes at a point if and only if A and its derivatives up to a certain order are zero 



[14], see for an explicit clear example the case of the massive scalar field below. We 
shall denote by A the set formed by the tensor A and its derivatives up to the required 
order such that, at any point x £ Vjv 

A\ x = 0^T(A)\ x = 0. (3) 

The set A essentially contains the derivatives of A appearing explicitly in the expression 
of T(A) for each case. 

The construction of T(A) generalizes the definitions of energy-momentum tensors 
of scalar and Maxwell fields |l^, Q and the traditional Bel and Bel-Robinson 
superenergy tensors of the gravitational field [12, 13, 14, 15, 16|, 17]. Many examples of 



T(A) will be given below. Note that the standard energy-momentum tensors of physical 
fields, even in the cases they are well defined, do not always have the dominant property 



[ 15| so they cannot be used as the tensor T{A) in general for our purposes. 

Let us try to see now under what conditions the physical field A will propagate 
causally. To that end, assume that S is an (N — l)-dimensional closed achronal set in 
Vjv and D + (S) its future Cauchy development (see |1C] for notation and definitions). 



As S is achronal, intD + (S) is globally hyperbolic and it can therefore be foliated 
by spacelike hypersurfaces = {t =const.}, where the gradient v = dt is timelike 
everywhere in int-D + (S). Take any point q £ D + (S), then K = J~{q) H D+(S) is 
compact and has boundary dK = S U H + (S), where S = S n J~(q) and H + (S) is the 
future Cauchy horizon of S. Using any of the superenergy tensors T(A) we define the 
superenergy integral 



w(t) 



[ T^- a -(A)v ai ...v amV = [([ T ai ~ am (A)v a2 . . . v^da ax | s , j dt' (4) 



Here J (Sj) denotes the causal past of St, rj the canonical volume iV-form on V/v ; 
dcr ai |s , the volume (N — l)-form on T,f pointing along v, and to the minimal value of 
t in K. 
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The procedure now follows the structure of a theorem for energy-momentum tensors 
in [ 10 1 . Differentiating, we get by Gauss' theorem 

dw 



dt 



w '(t)= T"-(A)v a2 ...v am da ai \j, t = (5) 



St 



= -/ _ T"™(A)v a2 ...v am da ai \ H+{ s ) - [T a i- a ™(A)v a2 ...v am da ai \s + 

JH+(S)nJ-(Xt) [ 1 Js 

+ I {[V ai T a ^(A)]v a2 ...v am + (6) 

JJ-{Y>t)r\K 

+T 0l - a -(A) [(V ai v a2 )v a3 ...v am +... + v a2 ... v am _, (V ai v a J] } V 

where da a \^ always denotes the outward-pointing volume (TV — l)-form on any (N — 1)- 
dimensional subset £. Given that T(A) has the dominant property, by (Q) and (|5|) we 
have 

w(t)>0, w'(t)>0. 

By (D we get that T a '- a ^(A)(V ai v a2 )v a3 . . . v am < M 2 T"™(A)v ai . . . v am on K for 
some constant M2, because V ai v a2 is bounded on the compact K. Similarily, we get 
constants M3, . . . , M m from the other terms so that putting M = M2 + . . . + M m and 
using the fact that the normal to H + (S) pointing outwards from K is future-pointing, 
we get by (||) that 

0<w'(t)<Mw(t)-f T a ^- a ™{A)v a2 . . . v am da ai \~ s +( [V ai T ai " a ™ (A)]v a2 ...v^rj. 

JS ' Jj-(Y. t )C\K 

If A = on S (we usually think of S as an initial slice), then due to the property (|3|) 
of the superenergy tensors T{A) we obtain 

0<w'(t)<Mw(t)+ [ [V ai T ai ~ am (A))v m ...v am r 1 . (7) 

Consider now the following divergence condition []: 

V ai T"-(A)v a2 . . . v am < hT a ^ a -(A)v ai ...v am (8) 



where h is a continuous function on D + (S). Then eq.(0) becomes simply 

< w'(t) < (M + H)w(t) 



1 This is enough for our purposes, but it can be obviously generalized to other more general conditions 
involving "powers" of the tensor T greater than one. On the other hand, (||) can be written in a form 
explicitly independent of the foliation, but this is not needed here. 
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where the constant H is the maximum of h on K. Integrating this relation (or using 
a rather simple GronwalFs inequality) and using w(to) = 0, one then arrives at 



w{t)=0. (9) 

This means that Tq___q(A) = on K (with e*o oc v) and hence T(A) = on K. For any 
point q £ D + (S) we can apply the above reasoning to the corresponding compact K 
to obtain T{A) = on K and hence we get that T(A) = on D + (S) by continuity, 
which is equivalent to A = on K due to (^). Note that changing future to past we 
also get A = in D~(S) and hence in ~D{S) = D+(S)UD~(S). 
Using the terminology of [[Tcfl , we have thus proven the following 

Conservation theorem: A\s = ==^ ^Ijj^sj = ^ ^ e superenergy tensor T(A) 
satisfies ^). 

This theorem is clearly related to the causal propagation of the field described by 
A, because if A ^ in any region arbitrarily close to D(S) then A will propagate 
in time according to the field equations but it will never be able to "enter" D(S), 
showing that A cannot travel faster than light. On the other hand, the universality of 



the superenergy construction |12|, |13|, |b|] is essential here. 

As an example in arbitrary dimension, consider a scalar field (f) with mass m sat- 
isfying the 'generalized Klein-Gordon' iV-dimensional equation (□ + m 2 ) cf) = k a V a <j), 
where k a is some continuous vector field. Define the tensor T{4>) as 

Tab = V a 0V fe - ^sw>V c (/>V c (/> + ^m 2 4> 2 g ab 

(this is the standard energy-momentum tensor if k a = 0, in which case T(<p) is 
divergence- free) . T{4>) is easily shown to have the dominant property in any dimen- 
sion and to be zero if and only if <f> = (iff V a </> = for the massless case) in open 



neighbourhoods [14|. However, at points x € V/v, we have [14| 



T(<f>)\ x = ^ <f>\ x = 0, V<f>\ x = ifm^O 
T((p)\ x = V0| X = O ifm = 

Thus, in this case the set A is given by A = {</>, V(/>} in the massive case and by 
A = in the massless one. In any case, T((ft) satisfies the divergence condi- 
tion (||) so that the theorem shows that <fi\s, V<ft\s = =^ 01 ^^ = if m / 0, 
(V4>\s = <t>\ D / S \ =const. if m = 0). In addition, due to the linearity of the 
generalized Klein-Gordon equation, this result implies that 4>\ e)(s) ^ s un iQ ue given the 
initial condition ^>|,5, V^l^. A similar proof works for the electromagnetic field using 
the energy-momentum tensor in dimension N. 
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Much more interesting is the case of the free gravitational field in ./V dimensions, 
because then no energy-momentum tensor exists (due to the Equivalence Principle 
and the non-localizability of the gravitational energy). However, one can use the su- 
perenergy tensors as follows. Let T be the generalized Bel-Robinson tensor given by 







Tabcd — C a e J Cb e df + C a e J Cb ec f — ~g a bC le J Ci ed f ~ -^9cdC a el ^ ' C be if + -g a b9cd,C hel1 Chelf , 

1~abcd 1~bacd 1~abdc 1~cdab 
T a acd = (N — 4:)T a cad , T a cad = — - \ CdefCd e ^ — -^9cdC hel ^ Chelf 



where C is the Weyl tensor. Then T has the dominant property WM (see [19, 1£ , [L 



for N=4, in which case T is completely symmetric and traceless) and vanishes at any 
point if and only if C is zero there [14]. Let g be any metric conformally related to g 
by 

g = e 2U g (10) 

Then C\ cd = C\ cd and 

V a C a bcd = V a C a bcd + (N- 3)U a C a bcd (11) 

which implies 

V a f a bcd = e~ 2U [V a T a bcd + {N- A)U a T a bcd + U a (T a cbd + T a dbc )- 

-(N - 4)U b T a cad - U c T a bad - U d T\ ac ] . (12) 



If V a C a bc d = 0, then VaT a bc d = and by fll2j ) we see that the divergence condition 
holds for T. This comprises all spacetimes conformally related to Einstein spaces 
{Rab = h-9ab)i including the important case of Ricci-flat metrics. Thus we have a 
simple proof of the causal propagation of free gravity (the Weyl tensor) in all such 
spacetimes, which are of big generality. The special case of General Relativity (where 
N = 4 and the Bel- Robinson tensor is divergence-free) was studied in fl6|] . 

In the case of General Relativity, for many fields one can construct a tensor T 
satisfying (||). Note that the energy-momentum tensors, which would be obvious can- 
didates, are not always well defined (as happens for massless fields of spin higher than 
one), and that in many cases, such as for instance spin-1/2 fields, the standard energy- 
momentum tensors do not satisfy the dominant energy condition |lq] . Therefore, we 
must construct some other tensor T in these cases in the same way as with the general- 
ized Bel-Robinson tensor above. As mentioned before, in fact one can always construct 
a superenergy tensor T(A) of any field A with the dominant property |l^, |l3|, O]. 
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Consider a massless spin n/2 field p>. This is described by a symmetric spinor field 
<-PAx—A n = <P(Ax-A„) (spinor indices are capital letters; see fl5| for notation) and the 
field equation is |15| 

V AlB '<PAx...A n =0. (13) 

For the neutrino (n = 1) and the photon (n = 2) this is a well-posed initial value 
problem, but for n > 2 there arise the Buchdahl algebraic constraints [jl8| involving 



the curvature for existence of solutions. We refer to [18, 15] for a discussion on these 
restrictions. Nonetheless, independetly of these restrictions, we can apply our method 
to prove causal propagation of the solutions in general. 
Define the tensor T(tp) by 

Tax-On = VA 1 ...A n ¥A' 1 ...A' n 

(with the usual identification a\ <-> between tensor and spinor indices (jDJ and 



where the bar denotes complex conjugation). As shown in [13], any sum of spinors 
squared in this sense has the dominant property. Note that for n = 2 this is the usual 
energy-momentum for Maxwell fields. By (|l3|), V ai T ai ... an = and as T(ip)\ x = <^ 
<p\ x = 0, the conservation theorem is valid for <p. The theorem also holds for more 
general field equations. In analogy with the generalized Klein-Gordon ^-dimensional 
equation, we can consider the simple generalization 



An 



for some continuous vector field k a . Then the divergence condition is also satisfied and 
therefore the theorem is applicable to such fields too. In any case, the uniqueness of (p 
on D(S) for a given initial condition <p\s holds. 

For massive spin n/2 fields, one can find consistent systems of equations (without 
any Buchdahl's constraints) for spinor fields (p, \ such as ]l~9[,[20| 



V Al B ^A 1 ...A n = VXA 2 ...A n B>, ^(Ax B 'xA 2 ...A n )B' = VtpA x ...A n (14) 

with /i, v and m 2 = 2[iv. This is a well posed initial value problem for all spins. 
Let 

T ai ...a n = a(pAx-An l PA' 1 A^-A' n + PXA^Ai-AnXAxA'z-A'n ( 15 ) 

with constants a, (3 > 0. Then T has the dominant property, T\ x = 44> <p\ x = = x\x 
and 

V ai T ai ... an = ((3V - ail)XA' 1 A 2 ...A n ^ A ' 1 A' 2 ...A' n + 

(3n(n — 1) a r>pf 
+ , + i\ e (A 2 *A 3 ( PA4...A n )DEFXA 1 A' 2 ...A' n + complex conjugates. 
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Here we have used eq.(Cl) of and denotes the Weyl spinor. For any a and f3, any 
product of two components of x an d V 3 with a bounded coefficient will be smaller than a 
constant times the time component T0...0 of T (so it is not necessary to make the natural 
choice j3u = a/x). Thus the divergence condition (||) is satisfied and (<p,x) propagate 
causally. The uniqueness of ((f,x) on D(S) given the initial condition (<p,x)\s follows 
again. 

It is important to remark that the fact that \i is a constant is fundamental in the 
above reasoning, because it allows to isolate the derivative of x ( e< l-(Cl) in which in 
turn is used to check the divergence condition. Thus, for instance, the case of a massless 
spin-3/2 field satisfies the Rarita-Schwinger equations which, due to the Buchdahl 
constraints mentioned before, must be modified in curved spacetime in an appropriate 
manner, see e.g.[|l], [21]]. The resulting equations in Ricci-flat spacetimes are formally 
similar to (|l4| ) with n = 3, but instead of the constant /i there appears the Weyl 
spinor contracted appropriately with x And this makes a big difference, for the 
derivative of x cannot be isolated now in general (it requires an invertibility condition 
on the Weyl spinor). Then, if we construct the corresponding superenergy tensor (15), 
which still satisfies the dominant and other required properties, it is impossible to 
check the divergence condition generically. This is in accordance with the well-known 
problems of the causal propagation of Rarita-Schwinger fields in curved backgrounds 
(see e.g.[p2}|), which are reminiscent of the classical Velo-Zwanziger causal paradox [23] 
appearing when coupled to electromagnetic and other external fields due again to the 
existence of some invertibility conditions, see e.g. [24]. 

To summarize, we have presented an extremely simple and fully geometric argu- 
ment for the causal propagation of physical fields in arbitrary dimensions. A sufficient 
condition, in many cases very easy to check explicitly, for causal propagation is given 
in terms of the divergence of an adequate superenergy tensor of the field. All our 
treatment is classical, although it may be also useful in quantum aspects of the re- 
ferred theories. In this sense, notice that if any field may be described by tachyonic 
particles then the divergence condition, which is related to the field equations, must be 
violated. In general, we can choose to use any T(A) with property (0) that satisfies the 
divergence condition if this is helpful, and in fact there usually exists an s-parameter 
family of such tensors for a natural number s [JIJ]. Nevertheless, from (|j) it is clear 
that only the symmetric part of T(A) matters and thus the divergence condition refers 
to this intrinsically unique superenergy tensor (see also the discussion in fl4|]). On the 
other hand, sometimes even the completely symmetric tensors will not be unique (for 
instance in the case of (ff^), where a and (3 remain arbitrary), but in these cases the 
divergence condition will hold for all of them and the conclusion about causal propaga- 
tion will of course be definitive. Furthermore, the condition (||) can be weakened, for 
it does not need to hold pointwise but only on integrated form (even weaker conditions 
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on the integral in ([?]) can actually be imposed for the conclusion to follow). One 
can also formulate the conservation theorem for sets more general than D(S), but we 
have restricted ourselves to a physically interesting situation. 

We believe that more sophisticated theorems for physical fields can be found using 
superenergy tensors to define norms and to do a priori estimates (as in p5[ ; with respect 
to this, see also the excellent and very recent review |26||). The fact that it is always 
possible to construct a superenergy tensor with the dominant property of any field is 
fundamental here. Closer connections with the theory of symmetric hyperbolic systems 
need to be investigated because it may well happen that our divergence condition 
(which is obviously connected to the field equations for ^4) provide some useful criteria 
for the hyperbolic symmetrizability of some field equations. All in all, we wish to stress 
the simplicity of our argument and the universality of its application. 

G.B. is grateful to The Swedish Natural Science Research Council (NFR) for finan- 
cial support. 
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